4j 


AM A A 


NASA-th-hj^j 


A 


,7 


AIAA-96-1454 

Static Test Induced Loads Verification 
Beyond Elastic Limit 

V. Verderaime and F. Harrington 
NASA Marshall Space Flight Center 
Marshall Space Flight Center, AL 


37th AIAA/ASME/ASCE/AHS/ASC 
Structures, Structural Dynamics 
and Materials Conference 
April 18-19, 1 996/Salt Lake City, UT 


For permission to copy or republish, contact the American Institute of Aeronautics and Astronautics 
1801 Alexander Bell Drive, Suite 500, Reston, VA 22091 



AIAA— 96— 1 454-CP 


STATIC TEST INDUCED LOADS VERIFICATION BEYOND ELASTIC LIMIT 

V. Verderaime,* and F. Harrington* 

NASA Marshall Space Flight Center Alabama 35812 


Abstr act 

Increasing demands for reliable and least-cost high- 
performance aerostructures are pressing design analyses, 
materials, and manufacturing processes to new and nar- 
rowly experienced performance and verification tech- 
nologies. This study assessed the adequacy of current 
experimental verification of the traditional binding 
ultimate safety factor 1 which covers rare events in 
which no statistical design data exist. Because large, 
high-performance structures are inherently very flexible, 
boundary rotations and deflections under externally 
applied loads approaching fracture may distort their 
transmission and unknowingly accept submarginal 
structures or prematurely fracturing reliable ones. A 
technique was developed, using measured strains from 
back-to-back surface mounted gauges, to analyze, define, 
and monitor induced moments and plane forces through 
progressive material changes from total-elastic to total- 
inelastic zones within the structural element cross sec- 
tion. Deviations from specified test loads are identified 
by the consecutively changing ratios of moment-to- 
axial load. 

Nomenclature 

/V = normal load, kips 
M - moment, kip-inches 
C = cross sectional limits, inches 
H = element thickness, inches 
w = element width, inches 
E = elastic modulus, ksi 
n = strain hardening exponent 
K = strength coefficient, ksi 
G = normal stress, ksi 
£ = normal strain 

Subscripts 
ty = tensile yield 
tu = tensile ultimate 
cy = compression yield 
N = normal 
M - bending 

1 = minimum measured strain 

2 = maximum measured strain 
k = zone number 
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L Introduction 

Increasing demands for more reliable and affordable 
access to space are promoting leaner and more innova- 
tive structural designs that invoke more reliance on 
experimental verification of their behavior and safety. 
This compelling shift raises concerns on how well veri- 
fication tests are implemented. Assemblies of large, 
high-performance aerostructures are inherently very flex- 
ible, and structural boundary rotations and deflections at 
externally applied loads approaching rupture may 
improperly transmit the binding verification loads and 
unknowingly reject a perfectly adequate design or accept 
a submarginal one. 

To sample this phenomenon, the slope and deflec- 
tion were calculated at the free-end boundary load on a 
hypothetical cantilevered beam (Fig. 1). Deflections and 
slopes are shown for the yield and ultimate strain limits 
calculated at the fixed end. The tangent of the free 
boundary slope 6 is a measure of the consecutively 
applied load decomposing from bending to bending-axial 
load ratio. Though the vertical scale is exaggerated, the 
slopes and displacements at the free end are propor- 
tional. 



Fig. 1. Boundary Load Deformation. 

The ultimate load to fracture was calculated to be 
twice the yield load, and the resulting ultimate strain at 
the fixed end was an order-of-magnitude larger than the 
yield. At ultimate loading, the predicted deflection was 
an 18“ slope resulting in over 30-percent bending-to- 
axial load ratio deviation. The adversity of this ratio to 
the verification criteria is dependent on how it feeds into 
and intensifies critically stressed regions and how it may 
change the failure mode. 

Though many codes and texts are available for pre- 
dicting inelastic strain responses from imposed inplane 
and bending loads, literature is mute on determining test 
combined loads from measured inelastic strains. A tech- 
nique was developed to analyze elastic-inelastic strains 
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measured from back-to-back surface mounted gauges to 
determine and verify transmitted test loads with specified 
applied loads. 

II. Elastic-Inelastic M aterials Model 

Modeling elastic-inelastic behavior could be very 
difficult 2 unless idealized into the simplest mathemati- 
cal expressions within the physical phenomena of the 
material and its application, such as the two parameter 
power expression/ 


<J= K e n , (1) 

where “n” is the strain-hardening exponent. In the linear 
elastic region, a < F^, the exponent is defined as n = 
1.0, and for <J> F^ , the strain-hardening exponent is 
calculated from uniaxial stress-strain data 

_ lo g(FjFty) £) 

lo gUjEty) 

The strength coefficient “K" is evaluated at the yield 
stress, which is the elastic-inelastic interface, 

k = £z = r ( 3 ) 

4 

These properties are directly applicable to normal 
stresses and strains without interpretation through 
theory. 


III. Structural Modeli ng 

A rectangular cross section element illustrated in 
Fig. 2 represents most structural components and 
regions as in beams, plates, and shells. Bending and 
inplane normal loadings are the most commonly 
measured components on this type of element using 
back-to-back strain gauges. They often may be suffi- 
cient to sample and verify the load transmission of more 
complex systems, including transverse shears. 

Surface mounted 



and strain, 

T N Tn 

eN [KwH] ' 

are uniformly distributed over the element cross section 
throughout the elastic and inelastic range. Because cross 
section planes are known to remain plane after elastic 
and inelastic bending, the inelastic bending strain also 
varies linearly along the thickness. However, the stress 
varies nonlinearly with Eq. (1) and the bending neutral 
axis is not expected to coincide with the cross section 
centroid. Since bending and axial strains are linear, they 
may be algebraically added as shown in Fig. 3(a). These 
combined strains are measured back-to-back at the ele- 
ment surfaces as £ 2 £/* where £ 2 is assumed to be 

greater than £j . Figure 3(b) illustrates the nonlinear 
bending stress distribution derived from the strain distri- 
bution using Eq. (1) and the shift of the bending axis to 
balance the moment. 



3(a) Combined Strain 3(b) Combined Stress 


Fig. 3. Strain and Stress Distributions Along 
the Element Cross Section. 


When t 2 < £fy, the combined strain distributions 
over the cross section are all in one elastic zone, and 
when £/ > , the distributions are all in one inelastic 

zone. The objective is to define the elastic-inelastic zone 
boundaries for all other measured strain combinations 
and to calculate their contributions to the total normal 
and bending loads. 

The net strain from any midplane y-distance along 
the element thickness in Fig. 3(a) is defined by the pro- 
portionality 


e y = y(0.5//+y)+£, , (6) 


The induced normal stress in Fig. 2, 


and its location is 


On = 


N 

wH ’ 


( 4 ) 


y = - (e y -e i)-0.5 H 

Y 


( 7 ) 
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The bending strain slope is 


-in + l 


7“ 


eiS) 

H 


( 8 ) 


The incremental normal load along the cross sec- 
tion thickness is the product of the induced stress and 
area. 


dN - W 0 y dy = wK(£ y ) n dy . 


Substituting Eq. (6) for the strain and integrating, all 
zone normal loads may be calculated from 

ic a 

(9) 

Cb 

where C a and Q, are the integration limits of a zone. A 
zone is bound along the y-axis by the surface measured 
strains, Ej and £% or by the material limit changes noted 
by E ty and e cy . Substituting the appropriate pair of 
boundary strains into Eq. (7). 

„ i , N a 

C a ,b = - (£a,b~£ l) - ~ , (10) 

provides the upper and lower integration limits of each 
zone. The yield strain may be tension or compression, 
where = -Efy is assumed for a symmetrical material. 
The normal load across the thickness is the sum of all 
the zone normal loads 


N k = 


wKy n 

*+l 


H ex 

— + — + y 
2 7 


"in + 1 


a r=ZNk- (ii) 

Bending strain along the thickness is given by e My 
= €y-£ft , and the neutral bending axis is defined by a 
zero bending strain = 0 ), from which €y = e N . 
Substituting into Eq. (7), the neutral bending axis is 

C M = - (£*-£,)- 0.5 H , (12) 

Y 

where the normal strain, e N , across the thickness is 
determined by substituting Eq. (11) into Eq. (5). Using 
Eq. (12), the incremental moment about the neutral axis 
is 


Mk = vvtfy 



£L 

7 


+ 7 


H £, 

H ex „ 1 


— + — +y 

~ + — + Cm 


2 y 

NJ 


n+2 

n+1 




_ 


(14) 


The moment about the cross section is the sum of all 
the zone moments. 


A#=XAf*. (15) 

A unit width, w = 1, is assumed for plates and shells 
from which normal loads and bending moments are 
defined by kips per inch and kip-inch per inch units, 
respectively. Using the strain distribution expression of 
Eq. (6), the stress distribution along each zone is given 
by 


Oy = K[ABS(e y )] n SGN(e y ). (16) 

Expressions in absolute form allow raising strains to 
odd powers. SGN ( ) is the signum function, which re- 
establishes the sign of the expression. If the function 
equals -1, the strain is negative. 

IY, — Normal Load and Moment Solutions 

As the induced normal and bending loads in Fig. 2 
increase, the strain distribution over the element cross 
section progresses from totally elastic to totally inelas- 
tic in four possible profiles. Given the values of the 
two measured strains, £, and e 2% the related profile is 
directly selected, and the zones and integration limits are 
decided as shown in Fig. 4. 

The induced combined normal and bending 
moment loads in each strain profile are resolved through 
a straightforward analytical routine summarized as fol- 
lows: 


• Using zone boundaries strains from selected Fig. 
4 profile, the integration limits C a b for each 
zone are defined from Eq. (10), and are sub- 
stituted into Eq. (9) to solve for the normal load 
Nix of each zone in the profile, 


dM ~ wo y {y-C M )dy - wK(e y ) n (y-C M )dy . (13) 

Substituting Eqs. (6) and (12) into Eq. (13) and integrat- 
ing, a zone moment about the neutral axis is calculated 
from 


• The sum of normal loads from all zones, Eq. 
(11), is substituted into Eq. (5) to obtain the 
profile normal strain. 

• The bending neutral axis C^is located using the 
total normal strain from Eq. (5) in Eq. (12) 
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(I) £ C y< £ l< £ 2< £ ty (II) £ cy <£ l <£ ty <£ 2 



(HI) e ! < £ cy <£ ty <£ 2 (IV) £ ty <£ l <£ 2 


Fig. 4. Strain Profiles Over Element Cross 
Section. 

• Eq uati on (14) determines the bending moment 
M lt in each zone about the neutral bending axis, 
which are summed in Eq. (15) to provide the 
desired profile moment. 

• Strain distribution £y and stress distribution a y 
are plotted over the thickness using Eqs. (6) and 
(16), respectively. 

This direct, though laborious, routine was reduced 
to a simple computer code requiring no detailed knowl- 
edge of its derivation. 

V. Normal-Bendin g Loads Program 

Profile (III), having the most zones, was solved 
and programmed as summarized above. Other profiles, 
having fewer zones, were adapted by resetting limits 
according to their zone boundary values and positions in 
the strain diagrams and applying them to their appropri- 
ate zones. 

•NORMAL/BENDING LOADS FROM STRAIN 
DATA 

•NMLFSD, Microsoft Quick Basic 


' MATERIAL PROPERTIES 

INPUT "ELASTIC MODULUS E=";ELM 

INPUT "YIELD STRESS Fty=";FTY 

INPUT "MAX STRESS Ftu=";FTU 

INPUT "STRAIN @ MAX STRESS Etu=";ETU 

ETY =FTY/ELM 

PRINT "TENSION YIELD STRAIN" ;ETY 
ECY=-ETY 

SHE=LOG(FTU/FTY)/LOG(ETU/ETY) 

PRINT "STRAIN HARDENING EXPO. n=";SHE 
K=FTY/(ETY A SHE) 

PRINT "STRENGTH COEF K=";K 
KO=K 

SHE0=SHE 

ECYO=ECY 

ETYO=ETY 

TEST DATA 

INPUT "ELEMENT THICKNESS H= ';H 
INPUT "ELEMENT WIDTH w=";W 

10 INPUT TEST MAX STRAIN E2=";E2 
INPUT "TEST MIN STRAIN E1=";E1 
IF E2<E1 THEN 

PRINT "MAX STRAIN < MIN STRAIN" 

GOTO 10 
END IF 

IF E2=E1 THEN El=0.975 
SLOP=(E2-El)/H 

PRO=3 

■USING PROFILE (III) (E 1 <ECY <ETY <E2) 

IF ECY<E1 AND E1<ETY AND ETY<E2 THEN 
ECY=E1 :PRO=2 

ELSEBF ETY<E1 AND E1<E2 THEN 
ECY=E1 :ETY=E1 :PRO=4 
ELSEIF E2<ETY AND ECY<E1 THEN 
K=ELM :SHE=1 :ECY=E1 :ETY=E2:PRO= 1 
END IF 

Nini=W*K*(E2 A (SHE+l)- 

ETY A (SHE+l))/(SLOP’(SHE+l)) 

ND32=W*ELM*((ETY A 2)-(ECY A 2))/(2*SLOP) 

NDI3=(ABS(ECY)) A (SHE+1HABS(E1)) A (SHE+1) 

NIII3=Nin3* W*K/(SLOP*(SHE+ 1 )) 

Nnrr=Nffli+Nffl2+Nm3 

PRINT "TOTAL AXIAL LOAD N=";NiilT 

sNm=Nnrr/w/H 

PRINT "AXIAL LOAD STRESS SN=";SNHI 

IF SNIII<FTY THEN 

ENm=SNm/ELM 

ELSE 

ENnHSNmxn 1/SHE) 

END IF 

PRINT "AXIAL LOAD STRAIN EN=";ENIII 
EMMffl=E2-ENm 
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PRINT "MAX BENDING STRAIN EM=”;EMMIII 

CMm= (ENm-E 1 )/SL0P-H/2 

PRINT "BENDING NEUTRAL AXIS CM='';CMIII 

Min 1 =-((E2 A (SHE+ 1 ))-(ETY A (SHE+ 1 )))/(SHE+ 1 ) 

Mm 1 =Mffl 1 *((E I +E2)/2+CMffl*SLOP) 

Mffl 1 =Mffl 1 +((E2 A (SHE+2))-(ETY A (SHE+2)))/(SHE+2) 
Mmi=Mmi*W*K/(SL0P A 2) 

MIII2=-((ETY A 2HE 1 A 2))*((El+E2)/2+CMm*SLOP)/2 
Mffl2=MIII2-K(ETY A 3MEl A 3))/3 
Mm2=MIII2*W*ELM/(SLOP A 2) 
MIII3=-((ABS(ECY)) A (SHE+1> 

(ABS(E1 )) A (SHE+ 1 MSHE+ 1 ) 

Mffl3=Mm3*((E 1 +E2)/2+CMm*SLOP) 
Mffl3=Mm3+((ABS(ECY)) A (SHE+2)- 
(ABS(El)) A (SHE+2))/(SHE+2) 
Mffl3=Mffl3*W*K/(SLOP A 2) 

MmT=Mm 1 +MIH2+Mm3 

PRINT "BENDING MOMENT M=";MHIT 

Rm=MmT/NmT 

PRINT'MOMENT/ AXIAL LOAD RATIO R=’\Rin 
■LIMITS 

CTY=(ETY -El)/SLOP-H/2 
CCY=(ECY-E1 )/SLOP-H/2 
ETY A=FTY/ELM 

' STRESS & STRAIN DISTRIBUTIONS 
OPEN "CLIP:" FOR OUTPUT AS #2 

PRINT "PROFILE=" ;PRO 

IF PRO=3 THEN 

YS=-.5*H: YF=CCY : MY=9 

M=MY-I 

DY=(YF-YS)/M 

EY3=0: SY3=0 

y=YS 

FOR 1=1 TO M 
EY3=(.5*H+y)*SLOP+E 1 
SY3=K*((ABS(EY3) A SHE))*SGN(EY3) 

WRITE #2,y ,E Y3,ENni,ETY A,S Y3 ,SNIII,FTY 

PRINT y,EY3,ENin,ETYA,SY3,SNIII,FTY 

y=YS+(I+l)*DY 

NEXT I 

END IF 

IF PRO=l OR PRO=2 OR PRO=3 THEN 

YS=CCY: YF=CTY: MY=9 

IF E2<ETY THEN YF=.5*H 

M=MY-1 

DY=(YF- Y S )/M 

EY2=0: SY2=0 


y=YS 

FOR 1=1 TOM 

EY2=(.5*H+y)*SLOP+El 

SY2=ELM*EY2 

WRITE #2,y,EY2,ENIII,ETYA,SY2,SNIII,FTY 

PRINT y,EY2,ENIII,ETYA > SY2,SNIII,FTY 

y=YS+(I+l)*DY 

NEXT I 

END IF 

IF PRO=2 OR PRO=3 OR PRO=4 THEN 

YS=CTY: YF=.5*H: MY=1 1 

M=MY-1 

DY=(YF-YS)/M 

EP1=0: SP1=0 

y=YS 

FOR 1=1 TOM 

EP 1 =(.5*H+y)*SLOP+E 1 

SP 1 =K*((ABS(EP 1 ) A SHE))*SGN(EP 1 ) 

WRITE #2 1 y,EPl,ENm,ETYA,SPl,SNHI,FTY 

PRINT y,EPl ,ENIII,ETYA,SP1 ,SNIII,FTY 

y=YS+(I+l)*DY 

NEXT I 

END IF 

CLOSE #2 
REM STOP 
CLS 

ETY=ETYO 

ECY=ECYO 

K=K0 

SHE=SHE0 
GOTO 10 

A sample printout of the program giving cross 
section characteristics derived from back-to-back strain 
gauge data is 

ELASTIC MODULUS E=? 10500 

YIELD STRESS Fty =? 38 

MAX STRESS Ftu =? 58 

STRAIN @ MAX STRESS Eyu =? .06 

TENSION YIELD STRAIN 3.619048E-03 

STRAIN HARDENING EXPO, n = .1505829 

STRENGTH COEF K= 88.59669 

ELEMENT THICKNESS H =? 1.4 

ELEMENT WIDTH w=? .74 

TEST MAX STRAIN E2 =? .02 

TEST MIN STRAIN El=? -.01 

TOTAL AXIAL LOAD N= 16.21604 

AXIAL LOAD STRESS SN= 15.65255 

AXIAL LOAD STRAIN EN= 1 .4907 19E-03 

MAX BENDING STRAIN EM= 1.850928E-02 

BENDING NEUTRAL AXIS CM= -.1637665 

BENDING MOMENT M= 17.29161 

MOMENT / AXIAL LOAD RATIO R= 1.066328 
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VI. Conclusions 

Experimental verification consists of two coher- 
ent, deterministic static test parts. Structural response 
within the elastic limit is verified with specified exter- 
nal loads representing maximum predicted operational 
environments. The ultimate factor of safety covers rare 
events, and its traditional and historical usage exerts the 
greatest influence on design and acceptance criteria. 

However, the order-of-magnitude larger strains, and 
therefore displacements, imposed by the ultimate factor 
of safety may distort the applied load transmission. The 
documented technique was developed to identify and 
assess verification load transfer discrepancy through 
back-to-back surface mounted strain gauge data, which 
is applicable throughout the elastic and inelastic range 
of the structural material. 


It is concerning that verification test results often 
report surface strain measurements to conform very well 
with predicted math models up to the yield point, but 
then unexpectedly deviate during the inelastic loading to 
premature fracture. Reasons offered are usually indefi- 
nite. Perhaps this suggested technique may extend the 
basis for a more definite test evaluation. 
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